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above the yrast line. The temporal dynamics of the process is investigated by means of a set of master equations 
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I. INTRODUCTION 

Recent experimental developments in laser physics promise 
to open the new field of laser-induced nuclear reactions in a 
domain of excitation energies that has not been explored so 
far. Efforts are under way |[T1 to generate a multi-MeV zep¬ 
tosecond pulsed laser beam at the Nuclear Physics Pillar of the 
Extreme Light Infrastructure (ELI) now under construction in 
Romania O and at the International Center on Zetta-Exawatt 
Science and Technology (IZEST) El. Eurthermore, theoret¬ 
ical proposals for the generation of coherent gamma-ray fre¬ 
quency combs at ELI have also been put forward 01. How 
will an intense laser pulse interact with a medium-weight or 
heavy target nucleus? The nucleus is a strongly bound system. 
Therefore, the laser-nucleus interaction is weak in compari¬ 
son to the laser-atom one. A reaction that differs significantly 
from the standard photon-induced nuclear reaction is expected 
to occur only if the photons in the laser pulse are coherent. 
Only then does the effective dipole width attain values in the 
MeV range, making it comparable with other characteristic 
nuclear energy scales. In this paper we accordingly consider 
the interaction of a strong coherent zeptosecond laser pulse 
with a medium-weight or heavy nucleus with mass number 
A. The pulse contains W = 10^ — 10^ coherent photons, the 
energy Ep per photon is several MeV, and the duration of 
the pulse is fi/a where a is of the order of several 10 keV so 
that hja 7^ 10“^^ s. Relevant questions then are: (i) How 
does the interaction of this laser pulse differ from laser-matter 
interaction in other areas of physics a? (ii) Which are the 
reactions that we expect to occur? The answers are obviously 
interesting in their own right and relevant for the layout of 
future experiments. 

In the present paper we provide a partial answer to these 
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questions by addressing the quasiadiabatic regime of the laser- 
nucleus interaction. In that regime, the process of photon 
absorption and that of nuclear relaxation are governed by simi¬ 
lar time scales. The paper follows the study of the perturbative 
regime by one of the authors lb] where the absorption process 
is much slower than nuclear relaxation. We hope to be able 
to address the sudden regime (characterized by the converse 
situation) in a future paper. A brief sunmiary of first qualitative 
results for the quasiadiabatic regime has been published in 
Ref. 0. 

Our approach is based on the master equation describing 
the excitation and relaxation of the nucleus under the influence 
of the external field provided by the laser. Multiple absorp¬ 
tion of coherent photons leads to nuclear excitation far above 
yrast. Setting up the master equation requires, therefore, the 
knowledge of the A-particle level density pA at high excitation 
energies and for large particle numbers A, expressed in terms 
of the single-particle level density pi . That is a challenging 
problem because near its maximum and for A > 100, pA is 
several tens of orders of magnitude larger than pi . An impor¬ 
tant preparatory step in our work has been the construction of 
a reliable approximation for pA in terms of pi ElISl. 

Use of the master equation renders possible the semi- 
quantitative study of the competition between photon absorp¬ 
tion, stimulated photon emission, photon-induced nucleon 
emission, neutron evaporation, and induced fission. In the 
absence of particle emission and fission, photon absorption 
would saturate at an excitation energy where the widths for 
absorption and for stimulated emission become equal. That 
is the case at the energy where the level density pA reaches 
its maximum. Neutron evaporation takes over at an energy 
below the saturation point. The combination of repeated neu¬ 
tron emission and continued dipole absorption by the daughter 
nuclei then produces proton-rich nuclei far from the valley of 
stability. Although the induced fission width is small in com¬ 
parison to all other widths, fission eventually terminates the 
reaction chain unless the laser pulse comes to an end before¬ 
hand. In the latter case, laser-nucleus interaction experiments 
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promise to shed light on nuclei at excitation energies far above 
yrast and far from stability. 

A qualitative description of the expected processes and a 
definition of the quasiadiabatic regime are given in Sec. 
The master equation and the transition rates are introduced 
in Sec. Ill This section also contains a semi-quantitative es¬ 
timate of the energies and time scales involved in the photon- 
absorption, neutron-evaporation and fission processes. Numer¬ 
ical results follow in Sec.lIVland the paper concludes with a 
discussion in Sec. El 


II. GENERAL CONSIDERATIONS 

Nuclei are bound by the strong interaction. As a conse¬ 
quence, even the interaction of the strong laser pulse defined 
in the Introduction with a nucleus is much less violent than 
the interaction of a medium-intensity optical laser pulse with 
an atom. We substantiate that st ateme nt with the help of the 
Keldysh parameter llTOll 7 = \/2mIuj/{eE) used in atomic 
physics. Here, m and e are mass and charge of the electron, 
respectively, / is the field-free ionization potential, and uj and 
£ are the frequency and electric field strength, respectively, of 
the laser pulse. The Keldysh parameter determines the dom¬ 
inant interaction mechanism in atoms. For 7 < 1 tunneling 
ionization dominates while for 7 1 the process is governed 

by multiphoton ionization. A small value of 7 corresponds in 
the optical regime to an electric field strength £ « 10^ eV/cm 
of the laser. Such a field distorts the Coulomb potential of 
the atom so strongly that electrons are set free. The nuclear 
equivalent of the Keldysh parameter is obtained by replacing 
m by the nucleon mass (m -> 2000 m), uo by the photon en¬ 
ergy (uj 10^ cj), and / by the binding energy of the last 
nucleon (/ ^ 10^ /). These substitutions increase the value 
of 7 by a factor 10^^. To return 7 to a value less than unity, 
the field strength would have to increase by that same factor 
10^^, i.e., £ would have to be of the order 10^^ eV/cm. This 
value roughly corresponds to the ratio of the binding energy 
of the last nucleon and the nuclear radius and, thus, to a dis¬ 
tortion of the nuclear potential that roughly corresponds to the 
above-mentioned distortion of the Coulomb potential in atoms. 
In comparison with a standard laser, the photon energy in our 
laser pulse is increased by six orders of magnitude. Such an in¬ 
crease falls short by a wide margin of the necessary increase of 
£ by eleven orders of magnitude. Therefore, the laser-nucleus 
interaction is governed by a value 7 ^ 1 of the Keldysh pa¬ 
rameter. That confirms our initial statement that even a strong 
laser pulse perturbs the nucleus only fairly weakly. In what 
follows we actually focus on sequential absorption of single 
photons within one pulse, i.e., on even weaker electric field 
strengths than required for simultaneous multiphoton absorp¬ 
tion. In contrast to atomic physics, in nuclei an additional time 
scale plays an important role, namely the nuclear relaxation 
time (see below). In the quasiadiabatic regime that time scale 
is comparable with the time scale for absorption of a single 
photon. Single-photon absorption and nuclear relaxation occur 
therefore more or less simultaneously. 

For photons in the MeV range, the product of photon wave 


number k and nuclear radius R obeys kR 1. Therefore, 
we consider only dipole processes even though quadrupole 
excitation is important for some nuclei at small excitation 
energies CD For ^ 1 coherent photons in the laser pulse, 
dipole excitation is governed by the effective dipole width 
A^Fdip. Here Fdip is the standard nuclear dipole width and is 
in the keV range. The amplification factor N applies in the 
semiclassical limit: For N coherent photons the matrix element 
for absorption of a single photon is proportional to ^/N, and the 
rate is proportional to N. After absorption of Nq photons the 
rate becomes proportional to (A^ — A'o)- Foi' N Nq wq may 
simply use N — Nq ^ N. For A^Fdip we use values around 
5 MeV. Coherence is vital in bringing the effective dipole 
width up to values that are comparable with other characteristic 
nuclear energy scales defined below. Without coherence, the 
probability for the processes investigated in this paper would 
be dramatically reduced. In the course of the reaction, up to 
A^o ~ b X 10^ photons may be absorbed. 

A further distinguishing feature of nuclei is that the nucleon- 
nucleon force is basically attractive, and that nuclei are self¬ 
bound. Nuclear properties are understood in terms of the shell- 
model potential (the mean field) plus the remaining “residual” 
nucleon-nucleon interaction. Because of the latter, distinct 
modes of nuclear excitation have the tendency to mix with the 
numerous other nuclear modes that are near the same excitation 
energy: the nucleus equilibrates. That property is absent in 
atoms. It qualitatively changes the treatment of the multistep 
photon absorption process. The time scale h/Tsp for equili¬ 
bration is expressed in terms of the spreading width Fgp, a 
manifestation of the residual interaction. The value of Fgp is 
known for low-lying modes with excitation energies of up to 
ten or twenty MeV where Fgp is of order Fgp « 5 MeV fV2\ . 

The ratio A^Fdip/Fgp relates the speed of dipole absorption 
with that of nuclear relaxation and, therefore, defines three 
regimes of the laser-nucleus interaction, (i) In the perturbative 
regime A^F^ip ^ Fgp, single excitation of the collective dipole 
resonance in nuclei dominates. That regime has been investi¬ 
gated in Ref. M- Consequences for future experiments were 
theoretically explored in Ref. ca. It was shown that if excited 
above neutron threshold, the time dependence of the nuclear 
decay is non-exponential, both in the neutron and in the gamma 
decay channels, (ii) In the sudden regime A^Fdip ^ Tgp the 
residual nucleon-nucleon interaction is irrelevant. Nucleons 
are excited independently of each other and are emitted from 
the common average shell-model potential. The potential read¬ 
justs after nucleon emission. If the duration time H/a of the 
laser pulse is sufficiently large, the nucleus evaporates. The 
sudden regime is so far unexplored, (iii) The quasiadiabatic 
regime A^Fdip ~ Fgp forms the topic of this paper. The nu¬ 
cleus (almost) attains statistical equilibrium between any two 
subsequent photon absorption processes. 

A quasiadiabatic process occurs when the energy El per 
photon is less than or comparable with the nucleon binding 
energy Eb of around 8 MeV. Then absorption of a single pho¬ 
ton does not lead to nucleon emission. Rather, the excitation 
energy is shared (almost) instantaneously with several or many 
other nucleons. The nucleus equilibrates. Another photon 
is absorbed, and the process repeats itself. Consecutive ab- 
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sorption of A^o ^ 1 photons leads to high excitation energies 
NqEl of the (almost) equilibrated compound nucleus. The 
excitation process terminates either with the laser pulse or 
when all nuclei have fissioned. For the duration time Hja of 
the laser pulse, significant emission of particles sets in only 
after several or even many photons have been absorbed. Pre¬ 
compound reactions show that the condition El < Eb is not 
absolutely necessary ifT^ . If equilibration is sufficiently fast, a 
quasiadiabatic process will occur also for El > Eb- 

With the nucleus close to equilibrium at all times, photon 
absorption leads to compound-nucleus formation at very high 
energies. That situation is distinctly different from Coulomb 
excitation imi or nuclear excitation by inelastic electron scat¬ 
tering Ga. In both these processes, modes of excitation very 
far from the equilibrated compound nucleus are formed at some 
excitation energy E. Such modes decay by particle emission 
and/or equilibration. If reached at all, the compound nucleus is 
formed at a much lower excitation energy than E. 

The equilibration mechanism being absent in atoms, our 
theoretical description differs from the strong-field approxima¬ 
tion in atomic physics CSIIISI and is related to the theory of 
precompound reactions ini. As mentioned in the Introduction, 
we describe the process in terms of a set of time-dependent 
master equations. In addition to dipole absorption, we take 
stimulated dipole emission, neutron evaporation, and fission 
into account. Induced particle emission is briefly addressed in 
SectioifVl 

We simplify the treatment in two respects, (i) We disre¬ 
gard spin altogether, (ii) We assume that after each photon 
absorption process, the nucleus attains full equilibrium. Both 
approximations simplify our treatment considerably. Without 
approximation (i), the number of master equations would be 
multiplied by the number of spin values considered. Without 
approximation (ii), the same would happen with regard to the 
number of configurations needed to describe equilibration at 
fixed excitation energy. While both approximations can easily 
be removed, the resulting increased complexity of the approach 
does not seem justified at this early stage of investigation and 
in view of the complete lack of experimental data. 

To justify approximation (i), we calculate analytically in 
Appendix [A| the distribution of nuclear spin values after Nq 
photons have been absorbed by dipole transitions by a nucleus 
with ground state spin zero. This is done for Nq ^ 1. The 
distribution peaks at spin J = and falls off very rapidly 
for larger values of J. Exact numerical results not based on the 
approximation A/q ^ 1 confirm this result, see Fig.[^ Even 
for a maximum number A/q = 5 x 10^ of absorbed photons, 
the nuclear spin does not significantly exceed 20 or so. This 
fact justifies our neglect of spin. 

Approximation (ii) implies that in the master equation we 
need not distinguish (as done in precompound reactions) the 
formation of n-particle n-hole states. Instead we work with 
the full set of equilibrated states at each excitation energy. To 
discuss the error made in approximation (ii), we observe that 
for a non-equilibrated system the mean number of n-particle 
n-hole pairs at fixed excitation energy is smaller than for the 
equilibrated system. That has two consequences. First, the 
number of states accessible for further dipole excitation is 



FIG. 1: (color online). Distribution of total spin values of the com¬ 
pound nucleus for several values of the number A'o of absorbed pho¬ 
tons. 


larger (because the exclusion principle blocks fewer states). 
Second, the mean excitation energy per particle or hole is 
larger, too. Therefore, neutron decay is more likely than in the 
equilibrated case (the number of available decay channels is 
increased). Both errors work in the same direction, leading 
within our approximation to an overestimate of the overall du¬ 
ration time of the excitation process. The resulting uncertainty 
can be compensated by varying the relative strength of photon 
excitation and neutron decay in the calculations. 


III. MASTER EQUATION 
A. Basic Approach 

With A the mass number of the target nucleus, we consider 
a chain of (n + 1) nuclei with mass numbers A — i where 
/ = 0,l,2,...,n, with an arbitrary cutoff at i = n. In nuclei 
with even mass number A, the states with spin zero at excitation 
energy between (k — 1/2 )El and {k + 1I2)El are grouped 
together and are jointly referred to as states (i^k). Here El 
is the photon energy and k = 1,2,.... The group of states 
with excitation energies in the interval 0 < E < (1/2)El is 
labeled (/, 0). The number of such states is determined by the 
level density p{A,E) of states with spin zero. For odd A we 
proceed analogously. For simplicity and in order to avoid the 
introduction of additional parameters we neglect the even-odd 
staggering of the ground-state energies as well as the spin- 
cutoff factor, and approximate the level density of spin 1/2 
states by interpolating between the values for the two neighbor¬ 
ing even A nuclei. In other words, we use the expression for 
p{A^E) valid for even A and given in Ref. ||9l indiscriminately 
for both even and odd A. We construct the time-dependent 
master equation for the total occupation probability P{i^kA) 
of these states as function of time t. The equation takes into ac¬ 
count dipole excitation by the coherent laser pulse, stimulated 
dipole emission by the same pulse, both for every nucleus in 
the chain, and neutron decay populating nucleus A — i — 1 af 
the expense of nucleus A — i. It is assumed 
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that within each group k of states in nucleus A — i, the 
occupation probability is equilibrated at all times and, thus, 
proportional to the total level density p{i, k) for that group. 
That assumption is characteristic of the quasiadiabatic regime. 
Neglecting the emission of charged particles we confine our¬ 
selves to a chain of nuclei with equal proton numbers. Likewise 
we do not take account of particle loss due to direct photon 
excitation of particles into continuum states. We address the 
ensuing limitations and possible corrections below. 

For the duration time 1/cr of the laser pulse, the states (i, k) 
are fed by coherent dipole excitation of the states {i^k — 1) and 


by stimulated dipole emission of the states (i, /c + 1), and they 
are depleted by dipole absorption exciting the states (i, /c + 1) 
and by stimulated dipole emission to the states {i^k — 1). Using 
Fermi’s golden rule we write the rates feeding the states (i, k) 
as k) with k' = k ±1. Here is the 

square of the transition matrix element. Neutron decay depletes 
the states (i, k) at the rate VN{i^k). Neutron decay of the states 
(i — 1, /c') in the nucleus with mass number ^4 -|- 1 — i feeds 
the states (i, k) with the rate T ]sf{i — l^k' ^ k). The master 
equation reads 


P{i,k,t) 


e(l/cr - k)[Wl^_^P{i, k-l,t) + k + l,t)]- P{i, k, t)[W^k-iP(.h k-l) + 

Wkk+iPi^k + 1)]| + y]rjv(i - l,fc' -)> k)P{i - - TN{i,k)P{i,k,t) . (1) 


Here, h = 1 and we have defined P(—1, /c, t) = 0. The dot 
denotes the time derivative, and 0 is the Heaviside function. 
The initial condition is P(i, /c, 0) = We require that 

neutron emission does not take place from the nucleus with 
mass number A — n and put r 7 v(n, /c) = 0 for all k so that 
nucleus {A — n) serves as a dump for the overall probability 
flow. For i = 0,1,..., n — 1 we have 

TN{i,k) = J2^N{i,k^k') . ( 2 ) 

k' 


Then Eq. (j^ implies ^ P(i, k, t) = 0, and the master equa¬ 
tion conserves total occupation probability. 

Induced fission is taken into account by introducing a diag¬ 
onal loss term —F/(i, k)P{i^ k^ t) in Eq. ([B. Here F f(i^ k) is 
the width for induced fission from state (i, k). Fission leads to 
a depletion of the total occupation probability. We do not keep 
track of the fission products. Therefore, fission eventually ter- 
nunates the laser-nucleus reaction. In the absence of induced 
fission and neutron decay (i.e., for r 7 v(L = Vf{i^k) = 0 
for all (i, k)), probability is conserved within the target nucleus, 
J]]^P(0,/c,t) = 0. For cr 0 (infinitely sharp laser energy 
with the laser pulse lasting forever) and TN{i^k) = 0 for all 
(i, /c), the target nucleus equilibrates, asymptotically {t oc) 
reaching the equilibrium distribution Peq(0, k) oc p(i, k) for 
all k values below and around the saturation energy. 


B. Transition Rates 

The transition rates have been defined, calculated, and dis¬ 
cussed in Ref. d). We present these rates here for the sake 
of completeness. Because of the approximations explained 
below in the calculation of the level density, particular values 
of the transition rates might differ from the ones presented in 


Ref. 131. Since h= Iwq use the expressions “width” and “rate” 
interchangeably. 


1. Dipole Transitions 


With Fdip the standard nuclear dipole absorption width and 
N the number of coherent photons in the laser pulse, the ef¬ 
fective dipole width for the ground state is given by A^Fdip. 
That expression holds for ^ 1 in the semiclassical approx¬ 
imation. With Fdip in the keV range and N ^ 10^ or 10^, 
the effective dipole width is of the order of several MeV. The 
value of NT dip serves as an input parameter for our calculation. 
Photon absorption of an equilibrated compound nucleus at ex¬ 
citation energy E is then governed by the effective absorption 
rate (A^r)eff(£^) = Nr^ip pi,cc{E)/pi,cc{Eg). Here /9acc(-B) 
is the density of accessible states and Eg is the energy of the 
ground state. The expression for (A^r)eff (F^) is valid as long 
as the number Nq of absorbed photons is small compared to N. 
That is the case for the calculations presented below. We equate 
(A^r)eff (the rate for population of the states (i, k) by 

dipole absorption of the states (i, k — 1)) with k) 

in Eq. Q. The stimulated dipole emission width for the inverse 
transition (i, /c) ^ {i,k — 1) is then given by detailed balance 
as {Nr)st{kEL) = NrQfi{kEL) p{i,k - l)/p{i,k). In this 
way, all dipole rates in Eq. <0 are determined. With p(i, E) 
the level density of nucleus {A — i) at excitation energy E, we 
approximate the density of states (i, k) as p(i, k) = p{i, kEp). 

The level densities p{A^E) are calculated using the expres¬ 
sions for the total level density of spin-zero states in nucleus 
^4 as a function of excitation energy E derived in Ref. @ as 
functions of the single-particle level density pi{s). A realistic 
linear or quadratic energy dependence of pi is considered. 


(1)/ X 2A (2). . 3A 2 

Pi P) = Pi P) = 


( 3 ) 




5 


Here V with 0 < 5 < V defines the range of the single-particle 
spectrum, and F is the Fermi energy. The single-particle en¬ 
ergies £i with i = 1, 2,... are obtained from the condition 
i = /q" pi(£^). The expressions ^ are approximately 

valid for A = 100 and A = 200, respectively. We have chosen 
y = 45 MeV and F = 37 MeV and keep these values con¬ 
stant throughout the neutron decay chain. These values also 
determine the chosen number of total bound single-particle 
states O, namely 148 for A = 100 and 360 for A = 200, 
respectively. 

When the number of nucleons is large, the method in fails 
to work at small excitation energies. In this region we use 
the Bethe formula 1^ . This is done for the first 65 MeV of 
the excitation energy for the case of medium-weight nuclei 
(A = 100) and the first 200 MeV for heavy nuclei (A = 200), 
i.e., for approximately 10% of the total relevant spectrum. 
For the density pacc(^) of accessible states we have used the 
Fermi-gas model described in Ref. ||9l and the same choices of 
pi as in Eq.J^. 

Figures [^nd give the widths for effective dipole ab¬ 
sorption, stimulated dipole emission, neutron emission and 
the induced fission rates for medium-weight (^4 = 100) and 
heavy (^4 = 200) nuclei, respectively. Calculation of the latter 
two rates is explained below. Figures and show that for 
both medium-weight (A = 100) and heavy nuclei (^4 = 200), 
the effective dipole absorption width (A^r)eff(F^) decreases 
slowly with increasing excitation energy E. Over a range of 
1000 MeV the decrease amounts typically to a factor of two. 
The stimulated dipole emission width {Nr)st{E), also shown 
in Figs. andstarts out at small excitation energy from a 
value much below that of (A^r)eff (F^) and increases monotoni- 
cally with E. It becomes equal to (A^r)eff (F^) at the maximum 
^max of the level density p{k). Significant dipole excitation 
above the energy F^max is not possible because stimulated emis¬ 
sion outweighs here absorption, {NT)st{E) > (A^r)eff(F^) 
for E > Fornax- For the two choices of pi in Eq. ^ we have 
Fornax = 533 MeV and F^max = 1200 MeV, respectively. 

We have tacitly assumed that the spreading width Fgp re¬ 
tains the value of « 5 MeV found at low excitation energy 
also at the much higher excitation energies relevant for our 
paper. As is the case for (A^r)eff (F^), the energy dependence 
of the spreading width is determined by the relevant density of 
accessible states. As shown for (A^r)eff(F^) in Figs. an dH 
that density changes only moderately with excitation energy. 
That fact validates our assumption. 


Here is the neutron binding energy of nucleus {A — i). 

In the second of Eqs. the lower bound is zero for E = 0, and 
k' is bounded by k'E l < kEL — BN{i). Eqs. Q are consistent 
with Eq. We simplify the calculations by assuming that 
^Ar(^) = Bn = 8 MeV for all i = 0,1, 2,..., n for the short 
nuclear chains considered. We thereby neglect the odd-even 
staggering of binding energies and level densities. These run 
in parallel and, therefore, largely compensate each other in the 
neutron decay widths. For both choices for the single-particle 
level density pi in Eq. (§, Figs. [2] and show that TN{k) 
rises steeply with excitation energy. While much smaller than 
(A^r)eff for small excitation energies, T^ik) becomes equal 
to (7Vr)eff 3tE = En and exceeds (A^r)eff for E > En> For 
the two choices of pi and T{E) = 1 we have En ~ 435 MeV 
and En ^ 1080 MeV, respectively. Both values are smaller 
than the corresponding values of F^max, F^max = 533 MeV and 
Fornax = 1200 MeV. We note that the crossing energies En are 
here larger than the ones presented in Ref. JTl, i.e., the neutron 
rates grow more slowly with excitation energy in the present 
case. That difference results from our treatment of the Fermi 
energy F. In Ref. 0, the Fermi energy changes with mass 
number A whereas in the present work, we consider F fixed. 
That is in accordance with the fixed value of B^ii) = Bn = 
V — F = S MeV. The difference F^max ~ Bn is sufficiently 
large in both cases and for both A = 100 and A = 200 to be 
physically significant, in spite of the inherent uncertainties of 
the calculation of level densities at high excitation energies. 
We conclude that neutron evaporation is the limiting factor in 
nuclear excitation by dipole absorption (provided that fission 
and proton decay can be neglected). 

In Eqs. 0 we have not taken account of the transmission 
coefficients T{E), i.e., the probability of formation of a resid¬ 
ual nucleus at excitation energy E under neutron emission. 
These obey 0 < T{E) < 1 and should multiply the integrands 
in Eqs. 0. Neutron decay is dominated by s-wave neutrons. 
Here and except for the slowest neutrons the transmission co¬ 
efficients are of order unity (231. We test the influence of that 
approximation (which overestimates the neutron widths) by 
multiplying all neutron widths by a common factor 1/2. The 
pairs of dark blue dashed lines in Figs. andshow the band 
of neutron widths for T{E) in the interval 1/2 < T{E) < 1. 
For T{E) = 1/2, En is closer to F^max- 

3. Fission 


2. Neutron Decay 

For the neutron decay rates we use the Weisskopf estimate 


According to the Bohr-Wheeler formula (211 modified by 
friction (22l . the maximum width for induced fission (reached 
at friction constant = 0) depends on excitation energy E 
essentially as 


f^N{h k) = 


1 


27rp(i, k) 


I 


(k+l/2)EL—BN('i) 


dE' p{iEhE') , 


r N{h k ^ E) = 


27rp( 


1 

k) J(k' 


{k' + l/2)EL 


{k'-ll2)EL 


dE' p{i El, E') . 

( 4 ) 


Tf{i,E) = {hwi/{2E))e^p{-Ef/T} . (5) 

Here, is the height of the fission barrier, uji is the frequency 
of the inverted harmonic oscillator that osculates the fission 
barrier at its maximum, and T is the nuclear temperature. With 
T~^ = {d/dE)lnp{i, E), the fission width increases very 
slowly with E. That is shown in Figs. [^and[^ where we present 
the calculated fission rates for huji = 4 MeV and 2 MeV and 
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for = 10 MeV (4 MeV), respectively. The two values of 
Ef correspond to ^4 = 100 (^4 = 200), respectively, the fission 
barrier being lower for heavy nuclei. We see that while perhaps 
competitive with (A^r)eff at low excitation energy and for very 
heavy nuclei, the fission width at higher excitation energy is 
never competitive with dipole absorption or neutron decay, see 
Ref. ||7I. We must keep in mind, however, that for a suffi¬ 
ciently long laser pulse the processes described by the master 
equation Q are ultimately terminated by fission. The induced 
fission width for E = 200 MeV for a medium-weight nucleus, 
for instance, corresponds to a half-life of approximately 5 zs. 
As shown in the following section by our numerical results, 
fission is therefore terminating the laser-nucleus reaction after 
several tens of zs. 

Actually, little is known about the dependence of cji and 
on ^ for excitation energies in the range of several 100 
MeV above yrast. We expect to decrease with temperature 
and, more significantly, as ever more neutrons are evaporated. 
That will speed up the fission process. We have not attempted 
to estimate the dependence of Ef on temperature and mass 
number. Such an estimate would require different techniques 
and is beyond the scope of this paper. 



E (MeV) 



E (MeV) 


FIG. 3: (color online). Same widths as in Fig.j^for A = 200 (360 
bound single-particle states). 


E while the effective dipole absorption width (A^r)eff (E’at) 
decreases slowly with E. The point En of intersection of the 
two curves depends on (A^r)eff(E^) only extremely weakly. 
For our estimate we therefore use the values En ~ 435 MeV 
and En ~ 1080 MeV given in Section [111^ for all values of 
(A^r)eff (E^) and T^iE) considered below. Since the induced 
fission rates are much smaller than the dipole absorption rate, 
we neglect fission in this first approximation. Furthermore, we 
approximate (A^r)eff (E^) by A^T^ip for all values of E. Be¬ 
yond the critical energy En, neutron evaporation dominates 
strongly, and it is practically impossible to excite nuclei to 
energies larger than En- Since En is smaller than the energy 
^max defined by the maximum of the level density, stimulated 
photon emission is neglected. 

Disregarding neutron evaporation we first determine the time 
and the number of photons needed to reach the energy in 
the target nucleus. We approximate the master equation for the 
target nucleus by the set of equations 


FIG. 2: (color online). Widths for effective dipole absorption (solid 
red line), stimulated dipole emission (long-dashed green line) (both 
for VFdip = 5 MeV), band for neutron emission (short-dashed dark 
blue lines) and band for induced fission (dashed-dotted light blue 
lines) versus excitation energy E for A = 100 (148 bound single¬ 
particle states). The two neutron decay widths are calculated from the 
Weisskopf estimate (T{E) = 1) and by taking for the transmission 
coefficients uniformly the value T{E) = 1/2. The two fission widths 
are obtained using the values huji = 2 and 4 MeV. 


C. Estimates 

Before presenting our numerical results we give some es¬ 
timates that show roughly what to expect. We estimate the 
time dependence of the mean excitation energies and the range 
of nuclei reached by the combination of multiple photon ab¬ 
sorption and neutron decay. We recall that the total neutron 
decay width T n{E) increases steeply with excitation energy 


p(o, k, t) = -7vrdipP(o, k, t) + 7vrdipP(o, k-i,t). (6) 

The initial condition is P(0, /c, 0) = ^/co- The solution 

P(0, k, t) = exp{-7Vrdipd (7) 

obeys ~ ^ times t. Considering 

P(0, /c, t) for fixed t as a function of k and using Stirling’s for¬ 
mula, we find that P(0, /c, t) has a maximum at /c = /Cmax = 
A^Tdipt with width \/k^^. The critical energy En with 
k^E^) = k^ = En/El is reached after absorption of 
Nq = kjsf E \/l^ photons. The time needed for the pro¬ 
cess is tjsf = Ejsf/{E lNY^i^). For A = 100, Ejy = 435 
MeV, El = 5 MeV, ATTdip = 5 MeV that gives No = 87 
photons and tjy = 12 x 10 “^^ s. The corresponding figures 
for A = 200 are Ejv = 1080 MeV, El = 5 MeV, NY^ip = 5 
MeV, A^o = 216 photons and = 8 x 10“^^ s. The laser 
pulse has the required length in time if a < 50 keV or 
cr < 22 keV, respectively. The speed of the process increases 
and the number of absorbed photons decreases as either the 
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photon energy El ov the dipole absorption width A^F^ip or 
both are increased. 

As the process described by Eq. ^ carries on, neutron decay 
actually depletes P(0, /c, t) and feeds P{l,k',t). Even though 
for energies below the neutron decay width Fat is much 
smaller than A^Fdip, the time needed to reach excitation energy 
En in the target nucleus is sufficiently large that nearly the 
entire occupation probability in the target nucleus is lost to 
neutron decay on the way. In the daughter nucleus, photon 
absorption is described by the same equation (|6]) save for the 
feeding term. Neutron decay is dominated by slow neutrons 
and, therefore, implies a loss of excitation energy of about 8 
MeV (the binding energy). Therefore, feeding of P(l, /c, t) 
by neutron decay occurs at k values that are about two units 
(assuming a photon energy of around 5 MeV) smaller than 
the ones for which the loss in P(0, /c, f) occurs. It takes about 
two absorbed photons to make up for that energy loss. The 
daughter nucleus decays in turn by neutron emission. The 
process repeats itself in the second daughter nucleus and so 
on. As a result, the maximum of the occupation probability 
moves to ever proton-richer nuclei and to ever higher excitation 
energies, eventually hovering near or at most one or two k units 
below Pat. 

An overestimate for the time needed to reach the nucleus 
with mass number {A — i) at energy Pat is obtained by disre¬ 
garding neutron decay in the target nucleus at energies below 
Pat, i.e., by using Eq. (|^ up to P = Pat, and by assuming that 
at P = Pat neutron decay takes over and populates the states 
(1, /cat — 2) in the nucleus {A — 1) (this scenario is consistent 
with El = ^ MeV). These are dipole-excited twice till the 
next neutron decay sets in, and so on. Since at P = Pat we 
have Fat (Pat) = PFdip(pAr), the part of the master equation 
describing the feeding of the first daughter nucleus by neutron 
decay of the target nucleus reads 

P(l ,/cat — 2, t) = —PFdipP(l ,/cat — 2, t) 

+ A^FdipP(0 ,/cat, t) . (8) 

The loss term accounts for photon absorption in nucleus (A—1). 
From here Eq. with i = 0 replaced by i = 1 takes over 
until the energy Pat is reached in nucleus {A — 1). That is the 
case after absorption of two photons. The process continues, 
alternating between Eq. ^ and the analogue of Eq. ^ for the 
second, third, ..., daughter nucleus. Except for the counting 
of k values, these combined equations all have the same form 
as Eq. ([^, and their solution, therefore, has the form of the 
right-hand side of Eq. (|7]). Hence 

P{i, kn, t) = exp{-Arrdipi} . (9) 

[kn + 3 ^)! 

The function P{i^knA) has its maximum at t = {k^ P 
3i)/(PFdip). In other words, while it takes 87 (216) pho¬ 
tons to reach Pat from the ground state of the target nucleus 
with mass number A = 100 (A = 200, respectively), it takes 
only 15 additional absorbed photons to move 5 mass units 
away from the line of stability. The additional time needed 
is 15/(PFdip). These figures show that once the threshold 
energy Pat for significant neutron evaporation is reached, the 


process quickly populates nuclei far from the valley of stability. 
The spread a of the laser pulse needed to reach the energy Pat 
is given by 

(Jo = PFdipPL/pAT . (10) 

Here Pat is independent of P^. With Pat = 435 MeV and 
Pat = 1080 MeV for mass numbers A = 100 and A = 200, 
respectively, Eq. ( [T0| ) defines a relation between the experimen¬ 
tal parameters cf, El, and PFdip. The formation of proton-rich 
nuclei sets in whenever cf < cfq. Clearly, these times are over¬ 
estimates because neutron decay is actually simultaneous with 
(and not subsequent to) photon absorption. After several ten 
zs, fission puts an end to the process, sooner for heavy nuclei 
than for medium-weight ones. 

IV. NUMERICAL RESULTS 

We calculate the time-dependent occupation probabilities 
P(i, k, t) for medium-weight (A = 100) and heavy (A = 200) 
target nuclei that interact with a short pulse of coherent MeV 
photons. We solve Eq. Q numerically for several choices of 
photon energy Pl, of the effective dipole width PFdip, and of 
the length (n-j-l) of the decay chain. Eq. Q is written in matrix 
form (including target and daughter index i = 0,1,..., n) 
as P = M.V with P{i,kA) P{k + ^(/i^max + 1)A) 

^max the maximum number of excitation steps considered. 
The matrix Af is independent of time and has both block- 
diagonal parts (fixed index i) describing dipole absorption, 
stimulated emission, neutron decay and fission, as well as non¬ 
diagonal feeding terms for the daughter nuclei [mass (A — i)] 
which are populated by neutron decay of their predecessors 
[mass (A — i + 1)]. The formal solution of Eq. 0 in vector 
form is V{t) = exp{Alt}P(t = 0). Because of the strongly 
varying rates and the large level densities involved (for instance, 
p(^max) — 10^^^ MeV“^ for A = 200), severe numerical 
problems may arise when one attempts to use standard linear 
algebra routines such as LAPACK ll^ for the calculation of 
the matrix exponential. 

Dedicated matrix exponential methods for solving systems 
of extremely stiff differential equations, i.e., equations for 
which the solving methods are numerically unstable, unless the 
step size is taken to be extremely small, have been developed in 
the last half-century ||25HT71 and applied, e.g., to nuclear fuel 
burnup calculations ll2^1 - [3Qll . In this work, we employ a new 
fast general-purpose semi-analytical matrix exponential solver, 
which is based on a combination of backward-stable matrix 
algorithms ED- First, the system of differential equations is 
brought to triangular form using Schur decomposition. Then, 
the triangular system is solved recursively using a combina¬ 
tion of Schur decomposition and back substitution algorithms. 
This yields the solution in analytical form as a combination 
of polynomials and exponential functions of time t. To verify 
numerical accuracy, we additionally carry out comparative cal¬ 
culations with other state-of-the-art solvers for the particularly 
demanding case of heavy nuclei with A = 200. 

Fig.shows the target occupation probability P(0, f) = 
P(0, kELp) in the absence of both neutron decay and fis- 
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FIG. 4: (color online). Occupation probabilities P(0, P, t) for El — 5 MeV, A = 100 and for dipole absorption and stimulated emission only, 
(a-c) Contour plots of the time-dependent occupation probability P(0, P, t) as a function of excitation energy E for (a) NT dip =1 MeV, (b) 
NT dip =5 MeV and (c) NT dip = 8 MeV. (d-f) Corresponding peak height (dashed blue line) (scale on the right y-axis), position (solid red 
line), and FWHM (long-dashed green line) of the occupation probability (scales on the left y-axis) for (d) PFdip = 1 MeV, (e) NTdip = 
5 MeV and (f) iVFdip = 8 MeV. 


sion for photon energy El = 5 MeV and for A = 100 and 
^rdip =1, 5, 8 MeV. Depending on the effective dipole width, 
the saturation energy Pmax = 533 MeV where stimulated 
emission limits photon absorption is reached after 50 zs for 
^rdip = 8 MeV (after 100 zs for A^Fdip = 5 MeV, respec¬ 
tively X The energy E^ above which Fat > (A^F)eff is reached 
for the case of A^Fdip = 5 MeV much later than the time es¬ 
timated in Sec. |IIIC| of = 12 zs. A neck-like artifact can 
be observed at the switching point P = 68 MeV in the calcu¬ 
lation method for the level densities p{i,k) and is best visible 
in Fig. 1^. Additional information is provided by value and 
position of the maximum and by the FWHM of the occupation 
probability shown in the lower part of the figure for the three 
cases considered. In accord with our estimates in Sec. IIII Cl 
the FWHM is proportional to ^/t and, after the first few zs, the 
peak height has a 1 /dependence. The linear dependence 
of the maximum on time holds only for the smallest of the 
three values of A^Fdip and does so only for small times. In all 
other cases stimulated dipole emission slows down the linear 
increase. The switching point at P = 68 MeV is also visible 
here in the shape of kinks in the three curves illustrating the 
peak position, height and FWHM of the occupation probability 
in Fig.|^. 


Fig. shows qualitatively similar results for a target nucleus 
with A = 200. The time dependence of the peak position, 
maximum peak value and FWHM again confirm our estimates 
in Sec. |IIIC[ except that the linear approximation for the peak 
position is only valid for short times t (except for the case 
^Fdip = 1 MeV). Reaching the saturation energy Pmax = 
1200 MeV requires a substantially larger number of absorbed 
photons than for A = 100. The switching point between the 
calculation methods for the level densities p{i, k) at E = 200 


MeV is also here visible, especially in Figs. & an d|gi. 

Laser excitation of heavy nuclei (A = 200) poses a numer¬ 
ically challenging problem, due to increased stiffness of the 
equation system. We have used the parameter set A = 200, 
Pa = 5 MeV, A^Fdip = 5 MeV, as testing ground for 
five state-of-the-art equation solvers that employ matrix ex¬ 
ponential methods: (1) the semi-analytical matrix exponen¬ 
tial method ED, (2) the Chebyshev Rational Approximation 
Method (CRAM) with partial fraction coefficients for approx¬ 
imation order 14 (291 [32l, (3) CRAM with partial fraction 
coefficients for approximation order 16 (29l[32l, (4) eigenvec¬ 
tor decomposition of matrix M. and (5) a modem scaling and 
squaring Taylor expansion algorithm (331. The five solvers 
were used to reproduce Fig.|^. The comparison shows that dif¬ 
ferently calculated P(0, P, t) values agree within an accuracy 
of 10“^. All methods produce some numerical artifact in the 
form of fringes on the upper side of P(0, /c, f) on an accuracy 
level of 10“^ — 10“^. We do not expect that such a level of 
accuracy can be attained in experiments in the foreseeable fu¬ 
ture and consider these fringes irrelevant. The semi-analytical 
algorithm used in the present work ED is the fastest one, yield¬ 
ing results in less than a minute, while the slowest scaling and 
squaring Taylor expansion method required more than two 
weeks for the same parameter set. 

While Fig.j^shows how the value of A^F^ip affects the speed 
of nuclear excitation. Fig. [^displays the influence of P^. We 
take El =1,5 and 10 MeV and use the same dipole absorp¬ 
tion width A^Fdip = 5 MeV throughout. The dependence of 
the occupation probability P(0, P, t) on photon energy Pa is 
shown in Fig. for the generic medium-weight target nucleus 
with A = 100 and in the absence of neutron decay. Depending 
on photon energy, the saturation region is reached after 400, 
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FIG. 5: (color online). Occupation probabilities P(0, P, t) for El — 5 MeV, A = 200 and for dipole absorption and stimulated emission 
only, (a-c) Contour plots of the time-dependent occupation probability P(0, P, t) as a function of excitation energy E for (a) PFdip = 1 MeV, 
(b) NT dip = 5 MeV and (c) PFdip = 8 MeV. (d-f) Corresponding peak height (dashed blue line) (scale on the right y-axis), position (solid 
red line), and FWHM (long-dashed green line) (scales on the left y-axis) of the occupation probability and for (d) NT dip = 1 MeV, (e) 
iVFdip = 5 MeV and (f) NT dip = 8 MeV. 


100 or 50 zs, respectively. A large photon energy speeds up the 
excitation process and may partially counteract the effect of a 
small dipole absorption width A^Fdip. For El = 1 MeV the 
excitation path is energetically more narrow than in the other 
cases. In summary, photon energy El and dipole absorption 
width A^Fdip, (i.e., the number of coherent photons in the laser 
pulse) jointly determine the time scale of the excitation process. 
We add few technical remarks: Due to the smaller excitation 
energy per step of the calculation, the case El = I MeV is 
more sensitive to where we switch from the Bethe formula to 
our method for the calculation of nuclear level densities. That 
point shows up as a small neck at ^ 68 MeV excitation energy 
in all three plots in Fig.|^ In addition, the numerical effort for 
El = I MeV is much larger than for higher values of El SiS it 
requires both a large matrix and, because of the slower exci¬ 
tation process, many more time points. The difference in the 
energy spacing also affects the contour plots, creating the false 
impression that the integral over occupation probability at any 
one time is smaller for small El . Numerically, for any time t 
the sum over all occupation probabilities P{0,k,t) equals 
unity with accuracy better than 10“^ for all three considered 
photon energies El- 

We now take account of neutron decay of the target and of 
three consecutive daughter nuclei, still neglecting fission. We 
solve the master equation for a chain of five nuclei with mass 
numbers ranging from Aio A — 4. We disregard neutron emis¬ 
sion by the last nucleus with mass number A — 4 which serves 
as a dump for the overall probability flow. The dimensions of 
the matrices M. are five times larger than for the parent nucleus 
only. In particular, for photon energy El = MeV the matrix 
dimension is 500 (1000) for an A = 100 (an A = 200) target 
nucleus, respectively. Fig. shows contour plots of the occupa¬ 


tion probabilities P(i, E’, t) for i = 0 (target) and i = 1, 2, 3,4 
(daughters) and for three parameter sets as indicated. In all 
cases, the final nucleus in the chain undergoes only dipole exci¬ 
tation with energies eventually reaching the saturation energy. 
Compared to pure dipole absorption, neutron emission is seen 
to broaden the distribution, cf. Figs. and[^. The occu¬ 
pation probabilities P(i, P, t) with i < 3 show that neutron 
emission comes into play early, slowing down the excitation 
process even at energies below E^, i.e., below « 435 MeV for 
A = 100 and 1080 MeV for A = 200, respectively. As soon 
as the neutron emission rates reach approximately 10^^ s“^, 
sufficiently many daughter nuclei are produced within the tens 
of zs time span of the laser pulse to strongly deplete the occu¬ 
pation probability of the target nucleus. A comparison of the 
probability distributions of A = 100 parent nuclei for photon 
energies El of 5 and 10 MeV energy [columns (i) and (ii) in 
Fig.[^ shows that higher excitation energies are reached (and 
correspondingly stronger neutron decay sets in) as the photon 
energy is increased. For t = 10 zs and for El = 10 MeV 
most of the nuclei have disintegrated to the dump {i = 4), 
while for Pl = 5 keV the target nucleus and the first daughter 
i = 0,1, 2 are still predominantly occupied at t < 10 zs. Com¬ 
paring columns (i) and (iii) of Fig. we note that for identical 
times t, higher excitation energies are reached in the heavier 
target. That is possible because of the much higher value of 
Pat = 1080 MeV for A = 200. As is the case for A = 100, 
neutron emission plays an important role long before the ex¬ 
citation energy E^ is reached. After 10 zs, the occupation 
probability of the target nucleus practically vanishes, while 
the occupation probabilities for nuclei with i = 2, 3 and even 
i = 4 are significant. 

Fission is taken into account by adding the fission width as a 
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FIG. 6: (color online). Contour plots of the time-dependent occupation probability P(0, P, t) as a function of excitation energy E for (a) 
El = I MeV, (h) El = 5 MeV and (c) El = 10 MeV. The results are for A^Tdip = 5 MeV, A = 100, and for dipole absorption and 
stimulated emission only. 


diagonal loss term in Eq. Q- We do not follow the fate of the 
fission products as their masses are much smaller than those of 
the target nucleus and of the first few daughter nuclei. For a 
one-to-one comparison, results are shown in Fig.j^for the same 
set of parameters as used in Fig ^ As expected from the decay 
rates shown in Sec. |IIIB[ Fig.j^hows that neutron emission 
is much faster than fission. For A = 100, fission is so much 
slower than neutron decay that the comparison of the results 
for the target (T) and the first three daughter (D1-D3) nuclei 
shows little difference between the cases without and with 
fission. Only for the last nucleus in the chain, where neutron 
decay is switched off, we do reach the time scale for which 
the fission rate produces significant loss. Since heavy nuclei 
(A = 200) have a larger fission width, some loss of probability 
can be observed already for the first daughter nucleus and for 
the daughters with i = 3 and i = 4 the occupation probabilities 
almost vanish. The effective loss of total occupation probability 
^ P(i, k^t) = 0 is displayed in Fig.j^for all three cases (i), 
(ii) and (iii). 


V. DISCUSSION 

The results of Section|IV]are obtained under neglect of the 
direct emission of nucleons by photoabsorption into the contin¬ 
uum. As shown in Ref. Q, such direct emission plays only a 
minor role for nuclei around A = 100 but is competitive with 
neutron decay for heavy nuclei. The effective charges of neu¬ 
trons and protons being nearly equal in magnitude, such direct 
photoionization is expected to produce neutrons and protons 
in about equal numbers. As a consequence, photoabsorption 
populates highly excited states not only in the chain of nuclei 
reached by neutron emission, but also in all nuclei with mass 
numbers {A — i) that lie between the valley of stability and 
nuclei in that chain. We have not attempted to calculate that 
process in detail. 

Fission ultimately terminates all the processes considered in 
this paper. Characteristic time scales are given in Fig.[^ The 
study of nuclei at high excitation energies and far off the line of 
stability is possible only when the laser pulse terminates before 
that characteristic time, i.e., if cr > Fy. It was mentioned 
above that we expect F/ to increase with increasing distance 


of the fissioning nuclei from the valley of stability. A reliable 
estimate would require a precise calculation of the height of 
the fission barrier Ef versus that distance. Termination by 
fission of the processes considered in this paper would allow 
to measure F f and, thus, to check such calculations. 


It was pointed out by A. Richter ll34l that at excitation ener¬ 
gies of several 100 MeV, the compound nucleus might undergo 
transitions in which excited states of the nucleon are popu¬ 
lated or in which subthreshold pion production occurs. How 
likely are such processes (which we have disregarded in this 
paper)? In the quasiadiabatic regime, the compound nucleus 
is near equilibrium at all times, and the answer follows from 
a statistical argument. The large number of equilibrium con¬ 
figurations at energy E must be compared with the very much 
smaller number of configurations where (almost) all the energy 
resides in a single mode (that of the excited nucleon or that 
of the pion). The situation is comparable to the emission of a 
fast neutron that carries almost all the excitation energy of the 
compound nucleus and leaves the residual nucleus in a state 
of low excitation energy. In comparison with slow-neutron 
emission, this process is suppressed by up to twenty orders 
of magnitude. We expect a large suppression factor also for 
the processes under discussion. However, each of these pro¬ 
cesses is easily distinguishable experimentally from any of the 
processes considered in our paper. In spite of their enormous 
scarcity, formation of any of the A resonances or subthreshold 
pion production might, therefore, still be observable. 


In conclusion, in the quasiadiabatic regime reactions induced 
by coherent laser light will cause the absorption of up to several 
100 photons. Such reactions offer the unique chance to explore 
the level density of the compound nucleus far above the yrast 
line, hitherto an unknown territory. The model we have used is 
plausible but not yet established. It is for experimental data to 
tell whether our present understanding is correct. The reaction 
will produce nuclei far from the valley of stability. Their decay 
processes will yield spectroscopic information not accessible 
so far. 
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FIG. 7: (color online). Contour plots of the time-dependent occupation probabilities P(z, E, t) with (from top to bottom) target nucleus {i — 0, 
label T) and four daughter nuclei (i = 1 to 4, labels D1 — DA) as functions of excitation energy E for A^Fdip = 5 MeV. Left column: Target 
nucleus A = 100, photon energy El = 5 MeV; middle column: Target nucleus A = 100, photon energy = 10 MeV; right column: Target 
nucleus A = 200, photon energy El =5 MeV. Please note the different color coding span for the case with target nucleus A = 200, which we 
chose for purpose of comparison with the next figure. 
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Appendix A: 

We calculate the distribution of spin values in the target nu¬ 
cleus after absorption of A^o photons. We take the direction of 
propagation of the photons as 2 ;-axis. In an unpolarized laser 
beam each dipole photon carries the ^-component angular mo¬ 
mentum ±1. Absorption of a single photon populates nuclear 
states with = ±1. These correspond to total spin J = 1. 
Absorption of two photons populates nuclear states with Jz = 
±2 (once each) and Jz = ^ (twice), corresponding to J = 2 


(once) and J = 0 (once). Absorption of Nq photons populates 
states with Jz = iA'o, T 2, zbA'o T 4,.... Inspection 
shows that the multiplicity of states with given \Jz \ is given by 
(^°) where 2fc = | J^| + ATo and = No/2, No/2 + l, ...,No 
for No even and k = {Nq + l)/2, {Nq + 3)/2,..., Aq for 
odd. We note that for Nq even (odd) only even (odd) values 
of Jz and, thus, of total spin J occur. The number Z{J) of 
such spin values is given by with 2k = J Nq. 

For A^o ^ 1 we use Stirling’s formula and approximate the 
difference by the negative derivative with respect to k. Thus, 


~ In ^^—4 exp{A^o In No - (1/2) {No + J) ln{No + J)/2- (1/2) {No - J) ln{No - J)/2} . 

k=iNo + J)/2 Ao + J 

I_^ 

By definition, Z{J) > 0. We note that Z{J) = 0 at J = 0. Z{J) is located at Jo = ^/Nq. We have used Jo <C Nq. For 

Inspection shows that Z{J) has a single maximum and drops A^o ^ 1 that condition is met by the solution ( |A1| ). 

off to relatively small values for J ^ Nq. The maximum Jo of 
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FIG. 8: (color online). Same as Fig.j^but including induced fission. 
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t(zs) 

FIG. 9: (color online) Loss of total occupation probability 
^^ P(z, k, t) due to fission for the three parameter sets used in 
Figs. El and[^ NT dip = 5 MeV in all cases, (i) target nucleus with 
A = 100, photon energy El = 5 MeV (solid red line), (ii) target 
nucleus with A = 100, photon energy El = 10 MeV (long-dashed 
green line), and (iii) target nucleus with A = 200, photon energy El 
= 5 MeV (short-dashed blue line). 




